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ABSTRACT 


The obiective of this studu is the enharrcement of existing experimental data 
on the aerodynamic characteristics of cylinders of elliptical cross section The 
model used in the stud^y has the fineness ratio of i,5 Tests were conducted at 
different angles of attack ranging from zero deai''-ee to rdn&ty degree in steps of 
ten degree and at three different Reynold Numbers in subsonic zone. The model 
was mounted horizontally in the 2-D wind tunnel The pressure distribution over 
the surface was determined with the help of a water muliiiube manometer and was 
used to calculate the lift, drag and moment coefficients for various angles of 
attack. The vortex shedding frequency was determind by a hot wire anemometer 
linked to dynamic signal analyser based on FFT technique 

Thie variation of drag coefficient, lift coefficient and moment coefficient of 
the model, with angle of attack is presented for three different wind speeds. 
Free stream turbulence level is also determiried for the operating range of speed. 
Similarly Strouhal no is plotted against the angle of attack for diffrent wind 
speeds. The calculated values of Strouhal nos. compare reasonably well with the 
reported data in the literature 
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CHAPTER 1 

IKTRODOCTION 

i.l INimiDUClTD^ TO ELUPJJC CmSS SECTION 

In recent yfears, elliptic cross sections have acquired a areat deal of 
significance, especially in the context of missile configurations and aircraft 
noses. These cross sections offer considerable aerodynamic advantages over the 
conventional circular ones, especially at high angles of attack. The actual cross 
section is not exactly elliptical as the upper surface is more convex than the 
lower half. However, since elliptical cross sections are the simplest geometrical 
shapes to approximate them, the present work has been done on elliptical cross 
section 

. The data existing in the published literatuf’-e on enquiry turned out to be meagre. 
While the country has embarked on an ambitious missile (I.G.M.D.P.) and aircraft 
<L.C.A. Sc L.T.A.) developement programmes, the importance of such experimental data 
can hardly be underestimated. 

Allen, H. Jullian and Perkins, Edward (195 i) have given a method for 

estimating the effects of viscosity on the forces and moments for inclined bodies. 
In their* recommended method, they have tried to correlate the cross flow in planes 
perpendicular to the inclined axis of the body to the flow around a cylinder and 
hence, the resulting cross forces are related to the section drag coefficient. The 
use of this method f^quires a knowledge of the section drag characteristics of 
the cylinder at different angles of attack. Jn the era of terrain mapping missiles, 
whose trajectory changes continously in response to topological changes, the angle 
of attack of cross flow will also vary continously, leading to a change in forces 
and moment acting on the missile. Since drag coefficient is different at different 
angles of attack, the mean pressure distribution will also be different. Similarly 
the vortex shedding frequency is a function of angle of attack arid hence the 
fluctuating pressure acting on the cylinder will also vary as angle of attack 
changes in response to changes in trajectory. 

The motivaiion for the study came with the successes notched up* by the 
country's space scientists in successfully firing missiles after missiles of 
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diverse capabilities A second motivating factor was the highlu sophisticated 
dvnanrac signal analyser recently acauir^ed by the aerodyriamics laboratory at 
cor iSider able expefise ar»d which can be errtplo*yed to measure the shedding 
freauency in turbulent wakes 

12 UTERATliRE REVEM 

Extensive research into literature reveals that very little data exists on 
cylinders of elliptical cross section. The field has attracted a few researchers 
but most of their work is of a perfunctionary naturG and fails to cover the 
subject exhaustively 

Noel K. Delany and Norman E. Sorenson^(iS53) have conducted detailed 
experiments to study the approximate variation of drag coefficient with Reynolds 
Numiber of several cylinders having different cross sectiorial shapes. They 
obtained data for cyliriders having cross sectional shapes of circle , ellipse of 
two fineness ratios (1:2 and 2:1), rectangle, diamond and two isosceles triangles. 
Data was obtair>ed for Reynolds Number as low as ilOOO and as high as 2.3 x IQ^. 
The vortex shiedding frequencies were measured and drag coefficients were 
determined. 

For the ellipse, they plotted a graph showing the variation of the drag 
coefficient and Sirouhal Number vs Reynolds number. For Reyrtolds Number less 
than 10" the Stouhal frequency was approximatly 0,2 . This value provides a 
confirmation for the results of the present experirrtental study 

V. J. Hodi and £. Wiiand^ <1970) have experimentally studied the aerodynamios 
of a set of two dimensional elliptical cylinaters with eccentricity of 0.8 and O.6.. 
during organised wake conditions. The data on Strouhal Nuriber, unsteady 
pressures and wake geometry is presented a.s a furiction of arigle of attack during 
static condition of the model. The effect of Reynolds Number on the fluctuating 
pressure is also examined. The results indicate dependence of the unsteady forces 
an Reunolds Number at zero angle of attack. The model has thirty pressure tap 
points distributed all along the surface. They have given graphs showing the 
variation of Strouhal Number with angle of attack. They have also plotted the 
variation of mean pressur^'e coefficients over its surface for four different 
angles of attack <« = . 30° 60® 90°). 
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^ 3 Scope of the Present Study 

In the present study the ob.jective is to measure the presssure distributian 
over an elliptic cylinder and to determine the lift coefficient, the draa 
coefficient, the moment coefficient at various wind speeds and varuina angles of 
attack The experiment was repeated to measure the Strouhal frequency again at 
different angles of attack and different wind speeds 



CHAPTER 2 


EXPERirCKTAL EQUPMEKT 

2.1 HIND TUNNEL 

The subsonic wind tunnel (known as 5 - D tunnel) at aerodynamics laboratory 
was used. It is a close circuit type having a rectangular cross section of i'x4^ 
with a length of The characteristics of the motor are as follows 



SCR 

DC 

Blower 

Drive 


Twin 


15 HP 


Input 

4Q0V. 

50 cycles 

3# 



57 - 98 - 57 

amperes 



output 

Armature 

il5V 

ilO ampere 

DC 


Field 

li5V 

5-7 ampere 

DC 


The motors are turned an by pressirtg the start buttons of respective units 
and their speeds are then controlled by manual speed ad.jlusting potentiometers. 
The maximum velocity obtainable is 42m/s. The wind tunnel has a plexiglass window 
with a circular groove in which a disk can r-otate The opposite side has a wooden 
trapdoor concentric with the circular gr*oove in the plexiglass. 

2:2 Pitot Static Tube 

For measuring the wind speed a pitot static tube was used. The pitot static- 
tube was placed far upstream in the flow • and the press»jr-e measurements were 
taken by the water manometer incliried at 30^ with the vertical. 
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2-3 FCX112 HicrrKnanom^ter 

The pitot static tube was connected to a FC0i2 micromanomeier It is 
supplied bu Furness Ltd. It can measure water height upto 19.9 mm and wind speed 
upto 64 m/s. The ‘values can be directly read off the digital display 

2-4 Mater Manometer 

It is an ideal instrument for pressure measurement. Inspite of its siirolicity 
it can measure pressures with great accuracy The multitube inclined water- 
manometer was used. The tubes are inclined at 30* to the vertical to increase the 
sensitivity. 

2^ Hot Mire Anetf^ometer : 

It has been used to measure the frequericu of vortex shedding by the 
elliptic cylinder. It uses a Hot-wire as a transducer and is designed to measure 
the instantaneous mass flow of gasses and liquid. It is ideal for measurircg high 
frequency flow fluctuations. It is based on the principle of convective heat loss in 
an electrically heated wire by the flow of gas or liquid surrounding the wire It 
basically determines the amount of power required to keep the temperature 
constant. The anemometer used is DISA 55A01 constant temperature type The hot 
wire is 10 % Rhodium and 90 % platinum It has a resistance of 40 fi/cm at 20* C. 
The length of the platirium wire is 1mm and diameter. 0 00762 mm , So the 1/d ratio 
becomes 132.23 which is greater than 100 as required. 

2 -6 Oscilloscxxae : 

It can measure frequency response upto 0 — 60 ko/s The signal from hot 
wire has been monitored coniinouslu on a .DBS 50 20A portable dual channel digital 
storage oscilloscope. It is incorporated with an 8 — bit CPU, with a maximum 
sampling rate of 1 MHz and effective storage frequency of 40 db Hz and a memory 
capacity of 1024 words • per channel. The saved waveforms can be displayed toeing 
fTtagnified upto 100 times . 
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2 7 Dynamic Signal Analyser : 

This is a powerful iristrument to analyse the incoming signals either in the 
frequency domain or time domain It works on the basis of FFT( Fast Fourier 
Transform). Such an analyser is useful for the analysis of the signals iri the 
range of a few millihertz to about a hc^idred kilohertz. An important property of 
FFT is that if the timte record has N equally spaced samples then the FFT 
transfof'mation from tiiT»e domain to frequency domain contains only N/2 equally 
spaced lines. This is because in the frequency domain, each line actually contains 
both amplitude and phase. The samples can foe averaged over a short time or long 
time. The obtained spectra can be stored in a floppy for later use or can be 
dir-ectly fed into a computer or a platter. Sensitive to small signals, it is operable 
in a dry environment upto 40^' C. In the present investigation a fP spectrum 
analyser is used. It can perform a few mathematical operations like FFT, inverse 
FFT, con.}ugate, separation into real and imaginary parts etc. Further, one can 
define one’s own mathematical functions which make use of either the present data 
or some stored data. 
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CHAPTER 3 

EXPERIMENTAL TEOHNIQyE AND MODEL MDUNnTPvIG 

3.1 MCTOEL 

The model was made of fiber reinforced plastic by means of casting. 
Initially a wooden model of elliptic cylinder was made. The major axis was kept at 
10.5 cm and minor axis at 7.0 cm. The length of the cylinder is 31 cm. With the help 
of the wooden model, the mould was made. This mould was used for making the FRP 
model. A core was used for making a through hole in the model. The model was then 

painted in order to get a smooth surface arvd reduce skin friction. In order to 

« 

measure the pressure, pressure tap points were required. Holes were made on the 
surface of the model by means of counterboring machine. The holes were made in a 
plane whichi was equidistant from either* of the end surfaces. The holes were made 
approximately at angular distances of iO^ . At the leading edge , holes were made 
at angular distances of 5^ in order to get higher accuracy. In all 44 pressure tap 
points were made. The position of the pressure tap points is given in fig A. 1 and 
table A.2( Appendix A ) . The pressure tap point at the leading edge is numbered 
one ar»d after that the points are numbered sequentially in an anticlockwise 
direction. Small metallic tubes of bore 1mm i^iere bent at 90^ and inserted into the 
holes after heating. These metallitic tiijes were inserted into the tap points after 
pushing them through the central hole in the model. These metallic tubes were then 
attached by means of rubber tubing to the manometer. All these rubber tubing were 
taken out through the central hole of the model. The model was fabricated at the 
Experimental Aerodynamics Division Laboratory.v NAL Bartglore and was kindly 
provided by Dr. M. Kr ishnamurthy^ Astt. Director* EAD, NAL Banglore 

3^ mUKTlHG OF THE MODEL 

The model was mounted on the wind tunnel with its axis torizontal and 
perpendicular to the flow. Two flat discs were attached to the model at either end. 
The discs were screwed irito the moctel so that they can be removed, when 
reQuired. One disc was made oF wood and fitted into a correspondirig blind hole 
macte into the wooden trapdoor of the wind tunnel. The other disc was made of 
plexiglass and was made to fit the circular groove in the plexiglass window. This 
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disc has a hole through which the rubber tubings connected to the manometer can 
be taken out Both the discs were made freely rotating in their grooves A wooden 
T was made to 'fit the hole in the plexiglass disc and was used for rotating the 
elliptic cylinder. Thet*^y the angle of attack could be changed. 

In order to determine the 0^ angle of attack, the wind tunnel was started at 
any speed. At zero degree angle of attack the stagnation pressure of pitot static 
tube and the presscr^e at the leading edge of the model should be exactly same. The 
position was determiried by rotating the model and checking the pressures till they 
are equal. A marker was attached to the plexiglass disc to correspond to zero 
degree and markings were made on the circular groove at angular distances of 10®. 

3-3 EXPERIhgNTAL TECHhgQUE 

Initially, the turbulence level in the wind tunnel was determined by means 
of a hot wire probe. The hot wire probe was connected to DISA type 55D26 signal 
conditioner & DISA type 55HiO CTA standard bridge from which thie signal was 
taken to a. DISA type 55D3i digital voltmeter and a DISA type 55D35 RtMS unit. The 
pitot static tube was connected to FC0i2 Furness micromanometer, Wiich directly 
gives the wind speed in thte tunnel. The initial mean value of the voltage of Hot 
Wire at zero wind was read. The temperature, cold resistance and the operative 
resistance were also read . Then the wind tunf*>el was started and the 
patentiometer was adjusted so that the speed in the wind tunnel is in thie same 
zone as the one in which experiments are to be conducted. At a fixed speed the 
mean and rms values of voltage of Hot Wire are read. Then the wind speed is 
changed and the process gets repeated. 

The model was mounted and the pressure taps were connected to the water 
manometer. The heights of the water columris in all the 44 tubes of manameter were 
read. Then the wind tunnel was started and the heights of wat45r columns in all the 
44 tubes of the manameter were recorded. The heights of the water columns in 
the U - tutoes connected to a pitot static tube were also recorded. Then the angle 
of attack was changed by 10® and the whole process was repeated. It was done till 
the angle of attack became SO®. Then thte wind speed was changed and the whole set 
of data was again collected. This was done for three different wind speed. 

For measuring the vortex shedding frequency hot wire anemometer was 
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used It was connected to dynamic signal analyser and an oscilloscope The hot 
wire probe was placed at a distance of 6a (a = semi-major axis of the elliptic body) 
downstream f rom the cylinder A vortex detectof*' was used to find the exact 
position of the Karman vortex street, so that the hot wire probe can be located 
there. After locating the probe in the vortex street, voltage fluctuation was 
studied on the oscilloscope If the probe is located correctly the voltage 
fluctuates wildly. The dynamic sigrial analyser is programed to take 15 sets of 
data. The vortex shedding freq/uency was also deter'^dned at three different speeds 
and different angles of attack (in steps of IG) 
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4.1 Data Reduction 


CHAPTER 4 

RESU,TS AND IMSCUSSIOt 


4,1^ Ttirbtilence Level 

The turbulence level in the irimd tunnel for different wind speeds was 
determined by means of the following formula 


Turbulence interisity 


firms „ A E 
1000 


®rm 5 =" value of voltage in millivolt 

E = mean value of Hot Wire voltage at a given i-jind tunnel setting. 

Eq = mean value of Hot Wire voltage at zero wind. 


4.1b Press(jre Coefficient 

The pressure at each tap point for each wind tunnel run and angle of attack 
setting, is measured by means of thte following formula ( as the height of the water 
column in the manometer is read in cm, all calculations are shown in C.G.S. units 
for greater accuracy and simple calculations^ 

P. = jc? g hj sini3 


where 

Pj = pressure at the ith tap point in dunes 
fl = density of water in gm/cc 
g =? acceleration due to gravity in cm/sec^ 
h. = < initial height of ith water column 

— present height of ith water column) 

0. = angle k#iich the water column tubes make with the horizontal 
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dynamic prassur-e head (Q) is measured by the fallowing for-rriyla 

Q = 9 hpg siniS 

rr dynarfiiG pressurs in dynas/crff'* 

Hps = difference between the heights of water column in the two 
U - tube connected to pitot static tube 


arms' of 


The stagnation pressure is the same as the pressure at the leading edge of 
the model when it is incTined to the wind direotian at ' angle of attack. 

The static prBSBUPB is determined by the fallowing forrfiula 

Ps = Pq - Q 

where 

Pg = static pressure of wind in dynss/cm^ 

Pq =5 stagnation pressure of wind in dynes/cm^ 

The stagnation pressure was taken to be same as the iDf^essure at the 
leading edge at 0"" angle of attack. The static p^rassura is determified by 
subtracting the dynaxrdc pressure from the stagnation pressure . Since the wind 
velocity was kept oonstant at all angles of attack the stagnati»on pressure was 
assumed to be constant at all angles of attack. 

The pressure coefficient at each pressure tap point is determined by the 

following formula. 

Pi -Pj 

_p. g 

wher« 

Cp. s= pressure cxseffioient at ith tap point 

4.io Drap 

The total drag force was determined by calculating the total force acting 
on the model along the direction of the wind speed. 

Let the major axis of the cross section correspond to X axis, positive 
direction being towards the trailing edge and the mino" axis correspond to Y axis, 
positive direction being counter clcokwise to x-axis. The scfuation of the ellipse is 
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determined by the following formula 



where 

X| = X co-ordinate of ith tap point 

y • = y co-ordinate of ith tap point 

a = semi-major axis of the ellipse 
b = semi-minor axis of the ellipse 

The slope of the tangent to the ellipse at the ith pressure tap point is 

r- - b 

rni = - ^ 

The slope of the normal at ith pressure tap point is given by 



where a/b is the fineness ratio. 

For this study the fineness ratio was kept at 1.5. so the above relation reduces to 
= i.5 

= tan”^<1.5 m^) 

where is the angle between the normal and the x-axis. 

The pressure at each pressure tap point acts along the normal at that 

point. 

Now to determine the force, a small element of the perimeter is taken. The 
force element acting at ith tap point is given by 

df 

This force elemeni acts on a surface element of unit length and width rA^l' 


where 
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= the distance of ith tap point from the centre of ellipse 
= (a^ 4- £ir}%p 

= the angular distance of line fining ith tap point from the -f-ve x 
axis, measur-ed -f-ve irt c-ounter clockwise direction. 


A«i' = 




i +1 


e* 


i -1 


Let a be the angle of attack. Then, the angular distance (0) between the 
dir^sction in which the force element acts and the wind speed direction is 
0 = 0^ — ^ 


Now this force element can be resolved into two components, along the wind 
speed direction and perpendicular to it. 

The component of the force element acting at ith point along the wind speed 
direction is df^^ cos<©^ - a) 
so drag force is given 

44 

D “ ^i GQs<^l - a.) 

i=l 

and the drag coefficient is 

c ■ — 1 ^' ■ 

■'D - Q 2a 


4 .id Lift CoeffiGierii 

The component of the force element acting at ith point perperidicular to the 
wind speed direction is df sinO^ - a) 
so lift force is given by 

' ■^44 

L = sirK^j^ - a) 

■ ■ i=l . ' ' 

and the lift coefficient is 
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^Sb Hcirr#jnt; Coefficierii 

The iTfonisrit about the oentsr of the allipse aoting on the cylinder due to the 
foro-a elerffeni dt^. acting at paint I- is given by 

, “4 “4 C 

dm^ = r- 3< uT^ 

= a df j_sinS^ cosS^-' — b df^c-osSj^ sinSj^' 


so the total moment acting on the cylinder 
44 

M = 2^ <P^ aS^') ( a sinSj^ cosS|' — boasS^ sinS|') in dyne-cm 

i=i 

The moment is taken as +ye if acting in clockwise direction. The moment 
coefficient is given by 


Cl 


M 


H 


Q <2a) <2b) 


4.if Sbrotihal hfejwfc^ 


The Strouhal number is given by = 


_ f d 


Uc 


where 

f 

d 

Ucio 


freauency of vortex shedding by the cylinder in Hz. 
2 b in meter 

free stream velocity in m/sec 


Hera for 'd', the minor axis has been taken in order to coiipare the results with 
the values obtained by Delany® 


4.1a ReyriOld Ho . 


Re No. 


a 


density of air in kg/m®(1.22S kg / m®) 


where 

fia 



free stream velocity in m/seo. 


?3 


d = 2 a in mater 

M = kinematic viscosity kg /m - sec(D.Q000i4S07 kg /rn - sec) 

45 Results . 

4 -2a TurtoulaiGe Level 

The turbulence lav'el in the operating zone turned out to be considerably 
high. The values are given in table a.i. It is more than 1.8% for all the three wind 
speeds, whereas ideally it should be less than 0.5%. .Since free strearf! turbulence 
is high, it will make the boundary layer unstable and induce transition from 
laminar- to turbulent boundary layer much earlier- than normal oircurfistances . 


4.2h Presstire Coef-Ficiartt 

The distribution of pressure coefficients over the surface of cylinder- far 
different angles of attack is plotted from fig. b.i to ta.iO . The three curves shcm 
the distribution of pr-essure coefficient for- different Reynolds no. If we study fig 
b.i, all the three curves appear to be quite symmetrically distributed about an 
axis passing throu^ x = 180”. This means that the distribution of pressure 
coefficients over- the top surface of the cylinder is quite similar to the 
distribution of pressure coefficients over the bottom sur-face of the cylinder. 
Since the pressure distribution is cwite identical over both the surfaces, the 
total lift force is nearly zero. The pressure coefficient remains constant from 
120® to 240®. This means that the flow separation htas taken place. Therefore the 
flow separates at about 120® from the leading edge for both the upper edge and 
the lower edge. The graph is quite similar to that obtained by Modi® in its general 
pattern. The Cp has a maximum value at S » 0*. In the case of Modi's graph ( fig 
6. a)® the pressure in the wake region is almost constant whereas in the present 
case it fluctuates sli^tly. 

If tsse study fig b.2 we find that a large suction pressure develops over the 
upper surface, between 0® to 100* whereas on the bottom surface the pressure 
remains positive only upto 40® and then becomes negative. However the value of the 
suction pressure on the bottom surface is very small. Therefore a very large lift 
force is generated. The pressure remains nearly constant from 120® to 240®. 
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ThsrsfQr-g the flow separ-aiion takes place at 120’ on the bottom surface and 60’ 
on the upper surface. 

Next we consider fig b.3. It is quite similar to fig b.2. Thera is a, very lar-ga 
v.alue of suction pressure on the top surface and small value of suction pressure 
on the bottom surface. However if we cOiT!p. 3 .re it with fig b.2,. we find that the 
suction pressure .acting on the upper surf .ace in fig b.3 is more than that in fig 
b.2 and therefore a larger lift force is generated at angle of attack -of 20’. The 
flow 5 ep.ar-.ates at about 120’ on the b-ottom surface and 130’ on the upper surface. 

Next is fig b.4. Here for Reynolds no. 2.29 x 10^ and 2.59 x 10^ , the 
suction pressure on the upper and the bottom surfaces are nearly equal and 
therefore they cancel out each other. There is a very small value of positive 
pressure acting on bottorri surface between 0 ’ to 60’ which is responsible for- 
generating senna lift. However the lift is substantially lower than in fig b.2 and fig 
b.3. The flow separates at about 120° on the bottom surface and about .30’ on the 
upper surface. These curves represent typical stall characteristics. However the 
curve for Re. no. 3.06 x 10^ differs and developes a large suction pressure on the 
top surface. This results in a large lift force. This behavior is difficult to 
explain theoretically. This may be due to experimental error. The graph an 
oomparision with Modi's ( fig 6 .b)® turns out to be quite similar. 

Fig b.3 corresponds to a = 40°. H^'e again flow separates on the upper 
surface at about 20° whereas on the bottom surf ace it separates at about 140°, a 
small value of positive prassur-e acting an the bottom surface gerierates a small 
lift. 

Fig b.S to fig b.iO are quite similar. !n all the cases the flow separates 
quite early on the upper surface and very late on the bottem surface. Far fig b. 6 , 
the flow separates at 140* on the bottom surface and 10* on the upper surface. 
For fig b.7 to fig b.lO the values are 160° on the bottom surface and 10* on the 
upper surface. The fig b.7 & fig b.lO show good resemblance to Modi's graphs <fig 
6.0 Se fig 6 .d respectively)®. 

in all tbiese curvej we also see that the maximun value of Cp occurs at the 
tap point ooresponding to the angle of attack. This is where the wind loses it 
entire kinetic energy and Cp attains it maximum value. 



4.2c Lift Coeff icisygj vs « 


Fig b.ii shoHs the variaiion of the lift coefficient vs angle of ^attack for 
the three Re. Nos. The three curves match perfectly but for one value at 30"" for 
Re. no. of 3.06 x 10^ . This anomaly beourries clear by viei^ing fig b,4 as 

observed, a large suoiion develops over the top surface resulting in a 

large lift force. This may be due to experimental error. 

The lift Gciefficient is nearly zero for zero degree angle of aitac.k. Ideally 
sfseakieig., it should be exactly zero. But because of the practical difficulty of 
iocatieig the exact position of the 0* angle of attack it ends up having small value. 

The lift Goefficierit increases rapidly for the next angle of attack 

arid continues to rise up to u - 20^4 After that the lift coeffioient 

drops suddenly representing typical stall charecterisiics Thereafter ■ the 
lift coefficient remains more or less constant , rising slightly to attain a 
local maxirnum at a = 50^. It finally becomBB zero for a = 90® , The curve 
for Re. no. 3. OS x 10 differs slightly by having a large value of, C|^ for 
a - 30® .This is due to large suctiorf pressure on the upper surface as 
evidenced by fig, b.4. 

4 2d larag Co&ffmimi^ vs m 

The graph in fig b.i2 shows the variation . The drag coefficient is 
slightly high for a as 0®This is because of flow separation as visible 
iri fi^ b.l. After this' the ct^ag coefficient drops and continues to have 
a small value for the riext angle of attack too .These two points refer 
to the zone where large lift force is generated .After ' that it' rises 

rapidly till a « 70®. After this point it rises more gradually and finally 

attains its roaxirrium value at « » 90®. The drag coefficients for zero angle of 
attack at all Re. nos. are approxiiiiaiely 0.4j whereas in the case of Delany^ they 
are more than 0.6. Delany^ has not given free stream turbulence level at which the 
experiment was conducted. The high turbulence level in the present study fffight be 
responsible early transition frcxii laminar to turbulent boundary layeri hence 
low value of drag caefficient. 
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4J2e Moment Coefficient vs a 


The 

relationship 

is 

platted 

in 

fig fo.l3 The moment coefficient 

is 

zero at 

zero degree 

of 

angle 

of 

attack because 

of symmetricity 

of 

pressure 

distributioTi 

over 

the two 

surfaces . As the 

angle of attack 

is 

increased 

,the morrient 

coefficient 

increases for <x = 

10* and reaches 

its 


peak value at oc = 20“^ .After that it drops rapidlu till 40^ .After that it 

remains nearly constant , attaining a local maximum at a = 60^ It also 
drops to zero for 90^ of angle of attack 

4.2f Strouhal Nuraber vs a. 

The spectral records for 0*^ and 10^ angle of attack for Re. no. 
2.75 X iO^ and 3.iS x iO^ do not give a clear peak. The curve remains more 
or less constant over a large range . The slight fluctuations might be due 
to voltage fluctuations and it is very difficult to determine the peak. 
The three spectral records are put at the end of Appendix C, the first 
one corresponding to Re. no. 2.75 x 10^ and it = 0* .The last two similarly 
for Re. no. 3.18 x 10^ and a = 0^ and a = 10^. The variation of the St. number 
vs angle of attack is plotted in fig.b.l4 .The value of Strouhal number has 
its maximum value at 0^ of angle of attack and after that it fTiare or 
less continues to fall and gives its lowest value at 90^. The results 

obtained more or less agree with the results obtained by Modi"'. The Sirouhal no. 

continue to fall as the angle of attack increases. However, the Re. no, 

dependence is not as marked as in the case of Hodi^. The Strouhal no. for zero 

degree angle of attack is approximately 02 as dotained by Delany . 
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5.1 Conclusions 

The pressure coefficient curves match quite well far all the Re. numbers. 
The pressure coefficient shows Re. no. independence. The pressure coefficient 
distribution over the surface is as expected. The graph for Cp vs & for the 
•angles of attack O'’’,. 30^'. £0'^ and 90^ .agree vary well with the results obtained by 
Modi"', Hodi's® graphs provide confirmatian for the correctness of method and allow 
us to corsclucfe that the results obtained for other .angles of attack are also 
reas-onably correct. The flow serapation takes place a.t all the Re. numbers. 

From the lift coefficient and drag coefficient curves we see that the 
cylinder is sub.iected to maximum lift force and minimum drag force when the arigle 
of attack is 20"^. The ideal angle of attack for generating large lift forces is 
between 10* to 20*. This fact is also corHPirmed by the rftoment coefficient curve. 
The maximum clockwise moment acting on the cylinder takes place at 20°. After this 
angle the cylinder exhibits typical .stall characteristics. The lift, drag and moment 
coefficients show Re. no. independence. 

The dynamic signal analyser turned out to be a ver-y useful instrument as it 
gives the frequency of vortex shedding directly without the botheration of 
counting number of peaks in the trace of the storage osoilloscc^je. Moreover, it 
averages the values over 15 reading, thereby increasing the acctr'acy. The 
Strouhal no. dea'-eases as angle of attack increases. 

5.2 Soogae fiy fuB"thar worfc 

One may study the variation of drag coefficient with Reynolds no. by 
oolleciing the data for large values of Reynolds Humber and can study the effect 
of tripping the boundary layer on the c^ag ooefficient. For this one may paste a 
small strip of sancte>aper near the leading edge of the model. 

One may also study the variatksn of the drag coefficient with changes in 
Strouhal number. 

The pressure drag coefficient is determined from the surface pressure data 
and the total drag coefficient is determined from wake measurement using a rake. 
This can be used to determine the skin friotiwi drag coefficient of the model. 
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Howev'er for this a very long wind tunnel will be necessary as in subsonic flow 
regime the downstream disturbarices effect the flow upstream. The one available 
presently in the aerodynamics lab is inadequate for this job. 

From the suT'Pbcg pressure data one can determine the separation and 
reattachrnent paints. This information can be confirmed by surface flow 
visualization using oil flow. 

In a nutshell, thfe field of elliptic cylinders is wide open to any student of 
experimental aerodynamics . 
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2m Model connected to water manometer 


3. Model connected to Furness micromanometer, 
dynamic signal analyser and plotter 
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TABLE a . 1 

Tanperature 301 K 
Cold Resistance 3 . 52 Ohms 
Operative Resistance 5.4 Ohms 


Wind speed 

E( V) 

Erms ( mV ) 

Turbulence 

( m/s > 



I nt ensi t y 

28 . 5 

6 . 97 


1.87% 

32.5 

7.09 

' 19.2 

1.88% 


38.4 


19.5 


1 .82 X 
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Re No. = 230000.00 



attack for an elliptical cylinder 



No. = 230000.00 







Fig b.l2 : Variation of drag coefficient with angle of 
attack for an elliptical cylinder 


Re No. = 2300C0.C0 
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appendix a 





TABLE AS 


Pressure 
Tap Point 

Angular 
Distancet deg) 

1 

0 . 00 

£ 

5.00 

3 

10.00 

4 

t S . 00 

s 

£0 . 00 

6 

£S. 00 

7 

31 . 00 

8 

36.00 

9 

, 42.00 

1 0 

48 . 00 

1 1 

S3. 00 

IS 

58.00 

13 

65.00 

14 

69 . 00 

IS 

75.00 

16 

83.00 

17 

90 . 00 

18 

98.00 

19 

105.00 

so 

111.00 

£1 

117.00 

£2 

128.00 


Pressure 
Tap Point 

Angular 
Oistance(deg ) 

23 

139.00 

24 

150.00 

25 

161.00 

£6 

171.00 

£7 

180.00 

28 

190.00 

29 

210.00 

30 

211.00 

31 

222.00 

32 

233.00 

33 

244.00 

34 

256.00 

35 

270.00 

36 

285.00 

37 

296.00 

38 

308.00 

39 

319.00 

40 

329.00 

41 

341.00 

42 

345.00 

43 

350.00 

44 

355.00 



APPENDIX B 



program progi input, output); 


const 

pi « 3.141S9E65; 

type 

vect « *rr*y Cl.. 443 of real; 

var 

1, prev, next,J : Integer; 

theta, alpha, b, e, c, d, s_9, s_tO,s_l! : real; 
theta_2, h : vect ; ~ 

table : array Cl.. 14, 1..443 of real; 

function sqi x real) -- real; 
begin 

sq : » X ♦ X ; 

end; 

function tan( x --real) -.real; 
begin 

tan -•» a inf xl/cosf x); 

end; 


begin 

for i 1 to 44 do 

readlni tabled, i3); 
read In ( alpha, c, d); 
b :« c; 
e •-« ( I . 9~c) ; 
for i :« 1 to 44 do 

readlnf tableCS,i3); 
for i t to 44 do 

readlnf tableC3,i3); 
for i :■ 1 to 44 do 

tableC4,i3 :* ( 180 . 0*arctan{ 1 .S*tan< {pi*tableE1 , i3 >/ 180.0))/pi) 
for i:*1 to 16 do 

tabIeCia,i3 • » tableC4,i3 +180.0; 
for 1 : *17 to £6 do 

tabled£,i3 tableC4,i3 +360.0; 

for i:»£7 to 3S do 

tabIeC1£,i3 -* tableC4, i3+360 . 0; 
for i:*36 to 44 do 

tabled£,i3 :* tableC4,i3 +540.0; 
for i :« 1 to 44 do 

tabIeCS,i3 ;« tableCIS, i3 -alpha; 
for i :« 1 to 44 do 

tablets, 13 :»{tableC3, 13- tablets, 13); 
for i t to 44 do 

tablet7, 13 (tablets, 13-e)/ b; 

for i :* £ to 43 do 


begin 

prev -• * 1 -1 ; 
next :* 1 +1; 

tablets, 13 (tabled, next3 - tabled, prev3)/E.0; 

end; 

tablets, 13 :« S.O; 
tablets, 443 :» 5.0; 



for i 1 to 44 do 
b#gi n 

theta :» { tabled, i3 ♦ pi)/ 180.0; 

tablets, i3 :«0.5* tableCB, i3 • «qrt( lO.BelO.S* 

sq( cosC theta)) + 49 * sq( sinC theta))) » pi/180.0; 

end; 

for i « 1 to 44 do 

tableC to, i3 -.■-tablets, i3* 

tablet 9, 13 * aln( (pi * tablet 5, 13)/180.0); 
for i .-*1 to 44 do 

tablet 11, i3 : --table t6 , i3 ♦ 

tablet 9, i3 * cost < pi * tablet 5, i3)/180.0); 
for i-. “t to 44 do 

tabletI3,i3 :■-<< table tS , i3 *table t 9,i3 * 
costt pi ♦ tablefIS, 13 )/180 . 0 )*{ain( { pi * table Cl , i3 )/1 80 . 0 ) ) *7 . 0 ) 
-(tablets, i3»tableC9, i3 ♦ sinCt pi * tableC 1 2, i3 )/1 30 . 0 ) ♦ 

(cos(( pi • table Cl , i3 )/1 80 . 0 ) )*1 0 . 5 ) ) ; 
for i:»1 to 44 do 

tabIeC14,i3 :«tabIeC1,i3 -180.0; 
s_9 : * 0 . 0 : 
s_10 :■ 0.0; 
s_l I := 0.0; 
for i :* 1 to 44 do 
begin 

s_9 : = s_9 + tablet 10, i3; 
s_l 0 := s_1 0 tablet 11, i3 ; 
s_1 1 := s_1 1 + tableC 13, i3; 

end; 

:* s_9/( 10.5* b) : 
s_l 0 :■ s_1 0/( 1 0.5*b) ; 
a__1 1 :■ s 1 1 / ( 1 0 . 5*b*7 . 0 ) ; 

writelnl ‘LIFT COEFFICIENT * ', s_9:15:5); 
writelnl 'DRAG COEFFICIENT « ', s_10:15:5); 
writelnt 'MOHENT COEFFICIENT = *. s_1 1:15:5); 

for i:*l to 44 do 

writelnt table CM, 13 •• 1 0 2 , table t? , i3 : 1 0 • E ) ; 


end . 
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Record at S6 . Im/s, lOdeg 
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Record at 2B . m/s, 50 deg 
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Record at 26. Im/s, BOdeg 






Record at 26m/s. 70 deg 
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Record at 26m/s, 90deg 







Record at O deg, 35.9m/s 
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